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We present experimental and numerical studies of localized terahertz surface waves on a 
subwavelength- thick metamaterial film consisting of in-plane split-ring resonators. A simple and 
intuitive model is derived that describes the propagation of surface waves as guided modes in a 
waveguide filled with a Lorentz-like medium. The effective medium model allows to deduce the dis- 
persion relation of the surface waves in excellent agreement with the numerical data obtained from 
3-D full-wave calculations. Both the accuracy of the analytical model and the numerical calculations 
are confirmed by spectroscopic terahertz time domain measurements. 
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Metamaterials have gained a great deal of interest in the 
last decade due to their important role as designer ma- 
terials with tailorable electric and magnetic properties. 
In this context, metamaterials provide a comprehensive 
tool box for the design of optical components with a 
very specific and in some cases even exotic electromag- 
netic behavior. From the conceptional point of view, 
the development of such optics requires the design of 
individual subwavelength-sized unit cells. When assem- 
bled properly, these subwavelength units stamp their 
electromagnetic characteristics on the effective macro- 
scopic properties of the resulting metamaterial. That 
way, it is possible to create optics with a very specific 
functionality. 

While free-space optics requires three-dimensional 
bulk metamaterials in most of the cases, it was shown 
that the surfaces of (thin) metamaterials (also termed 
meta-surfaces) can support the propagation of surface 
waves under well-defined conditions. Compared to the 
highly sophisticated methods that must be employed 
for the fabrication of bulk metamaterials, the tech- 
niques for the fabrication of most types of meta-surfaces 
are considerably less involved. As an intriguing anal- 
ogy, the observed surface waves on meta-surfaces are 
very similar in their properties to surface plasmon po- 
laritons on interfaces between a metal and a dielectric. 
In contrast to metals yet, the meta-surfaces can be 
specifically designed as to their effective material pa- 
rameters. In consequence, it is possible to tune the dis- 
persion characteristics and spatio-temporal behavior of 
surface waves on metasurfaces by design. 

For example, it has already been successfully demon- 
strated that the structuring of metal surfaces by small 
subwavelength slits or holes can lead to an improved 
confinement of surface waves. By this method, it was 
possible to localize surface plasmon polaritons near the 
metal surfaces at terahertz (THz) and microwave fre- 
quencies where metals almost behave like perfect elec- 
tric conductors P, 0, @] • The same concept was applied 
to the design of waveguides that consisted of structured 
metal sheets 3,|5|,|6|. Recently, in a similar approach, a 
waveguide structure based on complemetary split-ring 
resonators was proposed [7]. In an early experimen- 
tal work, it was shown that magnetic surface plasmons 
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Figure 1: (a) Schematic of one unit cell of the split-ring 
array, a — 41 pm, r — 17.25 pm, w — 7.5 pm, g — 3.25 pm. 
(b) Microscope image of a sample with grating structure. 



can exist on metamaterials with negative permeability 
[8|]. Moreover, it was discussed in a theoretical work 
that materials with simultaneously negative permittiv- 
ity and permeability can also sustain surface polaritons 
[9]. Facing the possibility of designing meta-surfaces 
with well-defined spatio-temporal and spectral proper- 
ties, it becomes obvious that such metamaterials offer 
the possibility to create compact, integrated plasmonic 
devices. However, to exploit the full potential of meta- 
surfaces with respect to the design process of tailored 
plasmons, it is necessary to understand the physical 
mechanisms in the interaction of meta-surfaces and sur- 
face waves. 

In this letter, we experimentally demonstrate the 
excitation of terahertz surface waves on thin mag- 
netic metamaterial films and provide a simple, intuitive 
model to describe the effective properties of the meta- 
surface. While previous experimental investigations in 
the THz regime focussed on structured metal sheets 

05 5 the magnetic metamaterials employed in this 
work consisted of a single layer of in-plane split-ring 
resonators (SRRs). The unit cell size of the structure 
was 41 pmx 41 pm, which corresponds to approx. A/7 at 
a frequency of 1 THz. The dimensions of the rings are 
shown in Fig.[TJa). The THz surface waves were excited 
by magnetic field coupling between the THz radiation 
and the magnetic resonances of the SRRs. It should 
be noted that the magnetic field of the guided mode 
is polarized normal to the surface in strict contrast to 
surface plasmons on metals. 

The meta-surfaces were fabricated by standard UV 
lithography. Details on the fabrication methods can 
be found in [10]. The copper SRRs were completely 
embedded in benzocyclobutene (BCB), a low- loss di- 
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electric with a refractive index of n = 1.63. The siUcon 
substrate which was used during the fabrication process 
was removed in a post-process. The resulting samples 
were free-standing metamaterial membranes with a to- 
tal thickness of approximately 50 pm. An exemplary 
microscope image of a sample is shown in Fig. [TJb). 
Since the dispersion relation of the localized surface 
waves is settled to the right of the light line of the 
surrounding material (air), incident THz waves can- 
not be directly coupled to the surface modes of the 
metamaterial structure. Therefore, we created an ar- 
tificial grating coupler by skipping every sixth row of 
split rings in x direction to convert THz waves to bound 
surface states. Due to diffraction at the grating struc- 
ture, the incident THz waves experience an additional 
momentum of kx = m x 27r/(6a) (m: order of diffrac- 
tion, 6a: grating period) parallel to the surface of the 
metamaterial. This allows that both the momentum 
and the energy are conserved in the conversion pro- 
cess if the frequency of the incident THz wave is be- 
low a specific threshold value given by the grating pe- 
riod 6a. In vacuum, for example, a parallel momentum 
kx = 27r/(6a) with a = 41pm corresponds to a fre- 
quency of ~ 1.22 THz. This means that below this 
frequency, the momentum of the first-order diffracted 
beam parallel to the surface is above threshold and the 
THz wave can couple to surface modes of the meta- 
surface. 

To gain a basic understanding of the propagation of 
localized surface modes on magnetic meta-surfaces, we 
developed a simple and intuitive model. For this pur- 
pose, we treated the metasurface as a thin slab waveg- 
uide with effective electric and magnetic material pa- 
rameters. In the model, the waveguide extended from 
x^y = — 00...00 and was confined in the z-direction 
between —d<z< -\-d such that the slab waveguide 
thickness was 2d. We assumed that the material was 
homogeneous and isotropic. To account for the reso- 
nant magnetic response of the split-ring resonators, we 
described the effective permeability /ieff using a Lorentz 
model 
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with resonance frequencies uoj , geometry factors Fj , and 
collision frequencies Tj [ll|. The effective permittiv- 
ity eeff was kept at a constant value that was equal to 
the permittivity of BCB. The material outside the slab 
was described by a permittivity eo and permeability 
jio (eo /io ^ 1 in air). To obtain the electromagnetic 
modes inside the waveguide, we considered a plane wave 



E = Eq exp [i{kxX + kzZ — cot)] 



(2) 



propagating in the x-z plane with a wave vector k = 
(kx^kz) and a frequency co. It should be noted that in 
this case kx is a complex quantity describing the phase 
propagation and absorption of the wave. The char- 
acteristics of the guided modes in the waveguide are 
mainly determined by the refiection at the boundaries 
between the slab material and the surrounding mate- 
rial. A stable mode pattern can only exist if the phase 
advance after two refiections at opposite boundaries of 
the material slab is an integer multiple of 27r. Since 
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Figure 2: (a) Dispersion relation of the localized wave, (b) 
Magnetic field normal to the metamaterial surface (z com- 
ponent) in the x-z plane and (c) in the x-y plane. 



a meta-surface only consists of a single layer of reso- 
nant elements, only mode distributions with a single 
maximum of the field amplitude in the center of the 
waveguide are meaningful solutions. These modes are 
obtained if the phase advance corresponds to 1 x 27r per 
round trip. Also, we only considered transverse elec- 
tric (TE) polarized modes (electric field vector points 
in y direction) to correctly account for the magnetic re- 
sponse of the SRRs in z-direction. For TE polarization, 
the Fresnel coefficient of reflection is given by 
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9okz ~\~ Meff^o2; 



(3) 



The wave number koz is determined by koz = 
[coiioiuj / cY — k'^Y^'^ and the sign of the square root 
is chosen such that |p| < 1. Finally, the condition for 
the dispersion relation of the localized modes can be 
written as 



arg(p^) ^Mkz = 27r 



(4) 



Note that p is a function of both co and kz. The solution 
of this equation is a dispersion relation kzioo)^ which can 
be converted to the more intuitive relation kx{oo) by use 
oikl^kl = eeff/ieff(^/c)^. The resulting function cj(/ca^) 
is indicated by the blue line in Fig. [^a). 

To check the validity of our simple model, we per- 
formed detailed numerical calculations using CST Mi- 
crowave Studio. Moreover, to obtain an accurate ref- 
erence to the experiments described later in this let- 
ter, we determined the dimensions of the fabricated 
metamaterial structure from a microscope image. The 
numerically calculated dispersion diagram in Fig. [2](a) 
(red crosses) shows two branches in the frequency range 
from to 2 THz in very good agreement with the re- 
sults derived from the waveguide model. It is notable 
that the waveguide model enables us to fully describe 
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Figure 3: (a) Measured and (b) calculated spectra of tera- 
hertz transmission through the sample. 



the subwavelength-thick meta-surface as an effective 
medium with effective material parameters. 

Besides the dispersion relation, we determined the 
magnetic field distribution of the guided modes from 
the numerical data. The results in Figs. [2jb) and (c) 
reveal that the magnetic field distributions of the two 
mode branches in the dispersion relation significantly 
differ near the plane of the split rings. At a given wave 
vector kx, the magnetic fields inside the rings and in 
the space between them oscillate in phase for the lower 
frequency mode, whereas they oscillate out of phase in 
the higher frequency mode. In the latter case, this leads 
to a stronger confinement of the fields in the immediate 
vicinity of the split rings since the fields cancel out with 
increasing distance from the plane. 

To support the results by experimental data, we mea- 
sured the transmission through the meta-surface us- 



ing a terahertz time-domain spectroscopy setup un- 
der normal incidence. Figure [3] displays the measured 
spectra together with the numerically calculated re- 
sults. The pronounced transmittance minimum at ap- 
prox. 1.1 THz for transverse magnetic (TM) polariza- 
tion (magnetic field vector in ^/-direction) results from 
the electric field coupling to the fundamental resonance 
of the split rings. For TE polarization (electric field 
vector in ^/-direction), one sharp resonance feature ap- 
pears at approx. 0.9 THz as well as two weak dips at 
1.15 THz and 1.2 THz. The frequencies of the transmis- 
sion minima are directly related to the corresponding 
eigenfrequencies of the excited surface waves with wave 
vectors equal to integer multiples of the reciprocal grat- 
ing vector, kx = mx27r/ (6a). The stronger minimum at 
0.9 THz corresponds to the excitation of a lower mode 
surface wave with m = 1 while the two weaker minima 
at 1.15 THz and 1.2 THz correspond to a lower mode 
surface wave with m = 2 and a higher mode surface 
wave with m = 1, respectively (see Fig. [2]). 

In conclusion, we have experimentally demonstrated 
the excitation of resonant terahertz surface modes on 
thin magnetic metamaterials by transmission measure- 
ments. The observed transmission minima could be 
related to the excitation of surface waves by direct com- 
parison to 3-D full- wave numerical calculations. In ad- 
dition, we developed a simple and intuitive waveguide 
model to describe subwavelength-thick metamaterial 
films as an effective medium. From this model, we could 
deduce the dispersion relation of the surface waves in 
excellent agreement with the numerical and experimen- 
tal results. The physical insight provided by this work 
is important with respect to the development of design 
tools for tailored surface waves on meta-surfaces. 
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